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We design a universal quantum homogenizer, which is a quantum machine that takes as an input 
a system qubit initially in the state p and a set of N reservoir qubits initially prepared in the same 
state 5. In the homogenizer the system qubit sequentially interacts with the reservoir qubits via the 
partial swap transformation. The homogenizer realizes, in the limit sense, the transformation such 
that at the output each qubit is in an arbitratily small neighbourhood of the state ^ irrespective of the 
initial states of the system and the reservoir qubits. This means that the system qubit undergoes 
an evolution that has a fixed point, which is the reservoir state ^. We also study approximate 
homogenization when the reservoir is composed of a finite set of identically prepared qubits. The 
homogenizer allows us to understand various aspects of the dynamics of open systems interacting 
with environments in non-equilibrium states. In particular, the reversibility vs or irreversibility of 
the dynamics of the open system is directly linked to specific (classical) information about the order 
in which the reservoir qubits interacted with the system qubit. This aspect of the homogenizer leads 
to a model of a quantum safe with a classical combination. We analyze in detail how entanglement 
between the reservoir and the system is created during the process of quantum homogenization. We 
show that the information about the initial state of the system qubit is stored in the entanglement 
between the homogenized qubits. 

PACS numbers: 03.65.Yz, 03.67.-a 



I. INTRODUCTION 



When a system interacts with a reservoir which is in 
thermal equihbrium then after some time the system is 
thermalized - it relaxes towards the thermal equihbrium. 
This imphes that the information about the original state 
of the system is (irreversibly) "lost" and its new state is 
determined exclusively by the parameters (temperature) 
of the reservoir. If the reservoir is composed of a large 
number, N, of physical objects of the same physical type 
as the system itself, then the thermalization process can 
be understood as homogenization: out of N objects (the 
reservoir) prepared in the same thermal state and a sin- 
gle system in an arbitrary state, we obtain N + 1 ob- 
jects in the same thermal state. This intuitive picture 
is based on certain assumptions about the interaction 
between the system and the reservoir, about the physi- 
cal nature of the reservoir itself and the concept of the 
thermal equilibrium. This picture is at the heart of the 
model of blackbody radiation, which triggered the birth 
of quantum theory in the seminal work of Planck. In 
addition, this same picture is very important in under- 
standing many processes in quantum physics as well as 
the fundamental concept of the irreversibility . 

In this paper we present a rigorous analysis of the 
above picture within the framework of quantum infor- 
mation theory. Specifically, we will consider a system, 
5, represented by a single qubit initially prepared in the 
unknown state g^g \ and a reservoir, R, composed of N 
qubits all prepared in the state ^, which is arbitrary but 
same for all qubits. We will enumerate the qubits of the 



reservoir and denote the state of the fc-th qubit as 

From the definition of the reservoir it follows that 
initially Cfc = ^ for all k, so the state of the reservoir is 
described by the density matrix 

Let U he a unitary operator representing the inter- 
action between a system qubit and one of the reservoir 
qubits. In addition, let us assume that at each time step 
the system qubit interacts with just a single qubit from 
the reservoir (see Fig. Moreover, the system qubit can 
interact with each of the reservoir qubits at most once. 
After the interaction with the 1-st reservoir qubit the sys- 
tem is changed according to the following rule (which is 
a completely-positive map) 



Tn 



(1.1) 
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FIG. 1. The scenario of homogenization with just three 
reservoir qubits involved. 

Let us repeat the interaction N times, that is, via a 
sequence of interactions the system qubit interacts with 
N reservoir qubits all prepared in the state The fi- 
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nal state of the system is then described by the density 
operator 



(N) 



(1.2) 



where Uk '■= U ® {®j^k ^i) describes the interaction 
between the fc-th qubit of the reservoir and the system 
qubit. This model of homogenization is very similar to 
the collision model since the system becomes homoge- 
nized via a sequence of individual interactions with the 
reservoir qubits. The interactions are assumed to be lo- 
calized in time (i.e., they act like ellastic collisions) 

Our aim is to in vest igate possible maps induced by 
the transformation (1.2) and describe the process of ho- 
mogenization. Homogenization means that due to the 
interaction U the states of the qubits in reservoir change 
only little while after N interactions the system's state 
become close to the initial state of the reservoir qubits. 
Formally, 



^N>Ns D{gf\0<5; 

Vfc,l<fc<iV D{^10<5, 



(1.3) 
(1.4) 



where £*(., .) denotes some distance (e.g., a trace norm) 
between the states, (5 > is a small parameter which is 
chosen a priori to the determine the degree of the homo- 
geneity and :— Trs[U g'g ® ^?7^] is the state of the 
fc-th reservoir qubit after the interaction with the system 
qubit. 
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FIG. 2. The 5 neighborhood of the reservoir's state ^ 
inside the Bloch sphere. After A'' interactions between the 
system and the reservoir the states of all reservoir qubits and 
the system qubit are contained within this (5-sphere. 



The conditions ( |l.3D and ( |l.4| ) can be represented using 
a geometrical picture. The Bloch sphere of unit radius is 



a representation of the state space of a spin-1/2 (qubit) 
system. The initial state p of the system qubit and the 
reservoir state ^ are represented by two (distinct) points 
of the Bloch sphere. We can image another sphere of 
the radius 5 centered at the point representing the reser- 
voir state ^ (in what follows we will call this sphere the 
(5-sphere). The task is to "shrink" the original Bloch 
sphere representing the (unknown) initial state space of 
the system qubit into the J-sphere. So we start with N 
reservoir qubits in the state ^ and the system qubit in 
an arbitrary state p and we end up with iV + 1 qubits 
within the (5-sphere centered at the point representing 
the original reservoir state ^ (see Fig. ^. 

We note that homogenization is closely related to ther- 
malization. There are however two main differences: in 
thermalization, (i) the state ^ of the reservoir qubits is 
not completely unknown, but is a thermal state, that is, 
a state diagonal in a given basis (interpreted as the ba- 
sis of the eigenstates of a one-qubit Hamiltonian) ; and 
(ii) the number of qubits in the reservoir is considered to 
be infinite for any practical purpose. Thermalization is 
studied in Ref. g]. 

Our paper is organized as follows: in Section || we 
show that quantum homogenization can be realize d w ith 
the help of a partial swap operation. In Section |lll| we 
show that the partial swap for qubits generates a con- 
tractive map on the system qubit with the fixed point 
being the initial state of the reservoir. This ensures the 
requi red conver gen ce of the homogenization process [see 
Eqs. (1.2) and (1.4)]. The uniqueness of the partial-swap 
operation is proved in Section IV . In Section |^ we es- 
timate the fidelity of the approximate homogenization 
map as a function of the number N of reservoir qubits 
and the parameter 5 (the precision of the homogeniza- 
tion), while in the Section VI we will analyze how the 
reservoir qubits become entangled as a consequence of 
thei r int eraction with the system qubit. In the final Sec- 
tion VII of the paper we address possible applications of 



the homogenization map. 



II. PARTIAL-SWAP OPERATION 

Let us start with the definition of the so-called swap 
operation S acting on the Hilbert space of two qubits 
which is given by relation H 



® \4>) = 10) ® IV). 

With this transformation 



(2.1) 



(2.2) 



after just a single interaction, the state of the system S is 
equal to the state ^ of the reservoir qubit; and the inter- 
acting qubit from the reservoir is lef t in the initial state of 
system. This means the condition (1.3) is fulfilled, while 
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the condition (1.4) is not — since recall that we want it 
to hold for all g^. 

In order to fulfill both conditions ( pT^ ) and (1.4) we 
have to find some unitary transformation which is "close" 
to the identity on the reservoir qubit, while it performs 
a partial swap operation, so that the system qubit at the 
output is closer to the reservoir state than before the in- 
teraction. The swap operator is Hermitian and therefore 
we can define the unitary partial swap operation 



P(r/) = cos 77 1 + i sin jj S 



(2.3) 



that serves our purposes. In what follows we denote 
sin rj = s and cos rj — c. 

In the process of homogenization, the system qubit in- 
teracts sequentially with one of the N qubits of the reser- 
voir through the transformation P{ri). After the each in- 
teraction, the system qubit becomes entangled with the 
qubit of the reservoir with which it interacted (for more 
details on the issue of entanglement see Sec. 0) The 
states of the system qubit and of the reservoir qubit are 
obtained by partial traces. Specifically, after the first in- 
teraction the system qubit is in the state described by 
the density operator 



(1) 2 (0) 



s'^ + tcs[tgfl (2.4) 
while the first reservoir qubit is now in the state 

^[^s'gf +c^^ + ^cs[gf,^]. (2.5) 

We can recursively apply the partial-swap transformation 
and after the interaction with the n-th reservoir qubit, we 
have 



An) 



Qs'=^'9s ' + ics Kofi's" ^^]' (2-6) 

as the expression for the density operator of the system 
qubit, while the 71-th reservoir qubit is in the state 



(2.7) 



Sinc e we a re in terested only in those terms in expressions 
(2.6) and (2/7) that are proportional to the operator ^ we 
can rewrite the above equations in the form 



n-l 

„(") _ „2 „2fc( 



gy^s-l^ c-^ + pill = (1 - c-)^ + p),Z , (2.8) 
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k=0 



and 



^2(n-l) 



+ 



(2.9) 



(2.10) 



Therefore, condition (L4) will be fulfilled for all k if and 
only if it is fulfilled for k = 1. This gives us a restric- 
tion on the parameter 77 that enters the partial swap; this 
restriction will be studied in Section 



III. HOMOGENIZATION IS A CONTRACTIVE 
MAP 

In this section we want to show that g^^'' ^ mono- 
tonically, for all para met ers 77 ^ 0. This means, in partic- 
ular, that condition (1.3) does not put any constraint on 
To show this convergence, we use the Banach theorem 
that concerns the fixed point of a contractive trans- 
formation. Let iS be a space with a distance function 
D{g,C)i then the transformation T is called contractive 
if it fulfills the inequality D{T[q\,T[£]) < kD{g,^) with 
< /c < 1 for all g, <^ S 5. A fixed point of the transfor- 
mation T is an element of S for which T[£] = ^. The Ba- 
nach theorem states that a contractive map has a unique 
fixed point and that the iteration of the map con- 
verges to it, i.e. r^[p] ^ for each g £ S. We note that 
contractive transformations within the context of quan- 
tum information processing have been recently discussed 
also in Ref. ||^. 

In our case S is the set of physical states, i.e. the set of 
all density matrices of a single qubit. The map, T, that 

we are considering is defined by g^g^ ^[Ps^] — ^?s^• 
We must show that the map is contractive, and that ^ is 
a fixed point of the map. 

We begin by finding the super-operator induced by 
the transformation U in the left-right form, i.e., as a 
linear operator acting qn_the space of trace-class oper- 
ators T(H) (see Ref. 



, cr,,, CTz 



We choose the operators 
(where cr„ are the Pauli matrices) as a ba- 
sis for T{TC), where H represents the Hilbert space of 
a qubit. In this case an arbitrary density operator of a 
qubit can be written as 



g^'-t 
^ 2 



(3.1) 



where |?7;| < 1/2. We can write a state that is an ele- 
ment of T(7i) in a vector form, i.e. g = (1, ly^;, ttJj,, ly^). 



Let ^ 



t ■ a = {l,tx:,ty,tz) be the state of the 



qubit in the reservoir. After the first interaction P with 
the first re serv oir qubit the system qubit evolves accord- 
ing to Eq.(2.6) with 77 = 1. This transformation can be 
described as 



In the next section, we are going to show that pill con- 
verges monotonically to the null operator as 77, —> 00. In 



this case, obviously g^^"^ ^ ^, so the condition (l.S) is 
fulfilled if the number of qubits N is large enough. In 
addition, as 77 increases, becomes more and more sim- 
ilar to since the commutator in (2.7) goes to zero; in 
other words. 



(0)i 



— 1 + (s^t + c^w) ■ a + ics[t ■ a^w ■ a] 



= -1 -I- [s^r + c^w- 2cs{tx w)] ■ S 



(3.2) 
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where we used the identity (7^(7; = Skit + i£jki<Jj, and 



(3.3) 



with j — X, y, z. Now we can express the transformation 

(0) (1) 
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(3.4) 

where T is the ma- 



or more formaUy, as Qg ' 
trix representing the super-operator acting on the hnear 
space T{Ti). If we express the matrix T as 



1 



0^ 
T 



(3.5) 



then it is easy to check that in our case Tt — c^t. This 
impHes that the state ^ is a fixed point of the map under 
consideration, i.e. = ^. The system state after n-th 
iteration then reads 



1 

-2^ 



n-l 

i=0 



3=0 



,2j 



(3.6) 



where for the last equahty we summed the geometric sum 
E"Jo(c^P = (l-c2")/(l-c2). Of course, c^" ^ unless 
c = cos 77 = 1. Numerically one can check that — > O, 
where O represents the zero operator. Thus gi^"'' — > ^. 
In what follows we prove this convergence for all values 
of the parameter rj. 

To prove that the map T is contractive, we must de- 
fine a distance function on S. Let us introduce the 
trace distance D{g,Lu) = Tr|g — uj\ and the vectors 
uj — (1, Vx, Vy, Vz) and f = w — v. For a qubit we have 

D{g,Lu)^TT\{w-v)-a\^TT\f-a\ ^ 2\f] (3.7) 

since the eigenvalues of the operator r • a are given by 
X± — ±|r|. In order to find the contraction parameter 
k for our transformation T we proceed as follows. From 
the Eqs. (3.5) and ( |3.7D we obtain 

D{T[g],T[uj])^Tv\f'-a\ = 2\f'\, 



where r ' — ' 

2 -> 
c r - 



-v' ^ sH+Tw-sH-Tv ^ T{w-v) = 
Tf = c^f - 2cstx f. Since 1^1^ < 1/4 and |r 'p = 
c^lfl^ -F4c2s2|rx f|2 = |r|2c2(c2 -|-4s2|i]2gjj^2^^^ where 
/3 < TT is the angle between the vectors t and r, we find 
that the contraction coefficient k — c. This last equality 
is due to the fact that \r'\ < \r\c. If c = cosr? < 1 then 
the map T is contractive and the convergence to the fixed 
point ^ is assured. 



IV. UNIQUENESS OF THE PARTIAL-SWAP 
OPERATION 



In what follows we will discuss the question of the 
choice of the unitary transformation, U, that describes 
the interaction between a system from the reservoir and 
the initial system undergoing the homogenization pro- 
cess. If both the system and the reservoir state are the 
same, the interaction should not affect either qubit, and 
this should be true no matter what the state of the sys- 
tem and reservoir qubit are. This implies that the unitary 
operator must satisfy the following two conditions: 



Tn (t/p®pC/t) = p, 
Trs{Up(E)pU^) = p, 



(4.1) 
(4.2) 



for any single-qubit state, p. Let us first discuss the case 
of pure states. If p represents a pure state then the con- 
dition (4.1) says that Up^pU"^ = p®£,i where ^1 needs to 
be determined. However from the second condition ( |4.2| ) 
it follows that Up ® pW = P where ^2 is unknown. 
Putting the last two results together we obtain that 



Up (g) pU^ = 



(4.3) 



for any p representing a pure state. From here it follows 
that the unitary transformation U acting on the joint 
Hilbert space Ti,^ —Ti^T-i must be of the form 



(4.4) 



where the parameter Lp is independent of the state IV)- 
Therefore, the action of the unitary transformation is 
fixed on the symmetric subspace of Ti,^ up to a phase 
factor e*'''. Neithe r of the two conditions (4.1) and (4.2) 



nor the condition (4.4) tell us anything about the action 
of the unitary transformation U on the antisymmetric 
subspace of Ti,^ . This means that action of U on the an- 
tisymmetric subspace is arbitrary. However, in the case 
of qubits the antisymmetric subspace is one dimensional, 
and we can proceed further. Because the antisymmet- 
ric subspace is one dimensional and invariant under the 
action of the unitary transformation J7, we have 



(3.8) [/(|V)|V-')-|V^)|V)) = e''(lV)lV^)-|V-')lV» , (4.5) 



where is a constant depending on U . Now the transfor- 
mation U is given by the equations (4.4) and (4.5) up to 
two constants </5 and 9. What we would now like to show 
is that these conditions require that t7 be a partial swap 
operator up to a global phase factor. This phase factor 
has no physical consequences. If we define the unitary 
operator U' to be 
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then equations (4^) and (|4.5| ) give us 



mm 



u[ pj:(e-¥^)/2 



1^)1^) 



Comparing these equations to Eq. ( |2.3| ), we see that C/' 
is just the partial-swap operator with rj — [ip — 9)/2. We 
can, therefore, conclude that in the case of qubits, the 
partial swap is the only possible o pera tor th at s atisfies 
the conditions of homogenization (11) and (12). The 
partial-swap uniquely determines yet another universal 
quantum machine the universal quantum homoge- 



V. APPROXIMATE HOMOGENIZATION 

In what follows we will analyze homogenization not as 
the limit of the infinite number of interactions, but as an 
approximate process after a finite number of steps. Let 
us s uppose that the parameter S from Eqs. (1.3) and 
(1.4) is fixed. This parameter characterizes our approx- 
imation. We will use the partial-swap evolution for the 
description of the homogenization. 

In the first step we gi ve a condition on the parame- 
ter T] of the partial swap (|2.3|). For our map T, we have 



,0 < 



D{gf\0- On the 
( |2.10D we know that D{Cn,^) < 



that D(ef \0 < 
other hand from Eq. 

D{Ci^_i,£,). As we have discussed earlier, we can adjust 
the parameter t] so that the condition D{£^[,£^) < S is ful- 
filled. Obviously, the distance D{^[,£^) depends on the 

initial state of the system, g^g \ and on r]. Therefore we 
have to determine the maximum value of ?7, for which 
the distance is less than or equal to S, independent (the 
universality condition) of the initial states of the system 
and reservoir. For a qubit the maximum value of trace 
distance is achieved for w = —t, corresponding to the 
situation in which the states are pure and mutually or- 
thogonal. The argument for this can be easily seen from 
a geometric representation of a qubit. In this case 



i?(ei,C) =2s2Tr|f. <j\^2s^ 



(5.1) 



since for a pure state |i| = i. From Eq. (5.1) we get the 
simple relation 



(5.2) 



sinrj < \/S/2 



The second step is to determine the minimum number 
of interactions, iV, that ensures for an arbitrary initial 
state of the system that the final state is in a sphere of 
radius 5 around the reservoir state ^. The worst case, i.e. 
when the number of necessary iterations is maximal, is 
intuitively the case when D(g^g\^) is maximal. In sec- 
tion III we proved the convergence of the system state 



to ^ for any rj 0. Therefore we are sure that such an 
TV exists. As was just discussed in the previous para- 
graph, the distance D{g^g\^) is maximal, when the two 
states are pure and mutually orthogonal. Moreover, our 
transformation T doesn't change the commutation rela- 
tion, which is initially equal to zero, i.e. [g^^\ £] — 0, for 



all N. Introducing w 
obtain 



-t for the commuting states we 



-l + (l-2c2^)f-a. 



and for the distance we find 



t) • ct| = 2c2^Tr|r- a 



(5.3) 



(5.4) 



This distance is maximal if we fix N and maximize over 
all g^'p and ^. Again, since |t| = ^ for pure states, we 



obtain the distance D{g\!^\£,) = 2c^'^ = 2(cosr;)2^. If 
the parameters and s in the experession (p^) are such 
that sin rj = \/5/2 then we can find the lower bound Ns 
on the number of reservoir qubits which are necessary to 
achieve the homogenization with a required fidelity 



N >Ns 



In (5/2 



ln(l - 6/2) 



(5.5) 



Both bounds on the parameters rj and N are completely 
determined by the parameter 5. After performing N it- 
erations, + 1 qubits are in states belonging to the 5 
neighborhood of the initial state of the reservoir, no mat- 
ter what the states f and g^g^ were. 

We see that if we fix the number of reservoir's qubits 
N , then t he o ther two parameters are determined by the 
relations (5^) and (p^). 



VI. ENTANGLEMENT VIA HOMOGENIZATION 

In spite all the progress in the understanding of the na- 
ture of quantum entanglement there are still open ques- 
tions that have to be answered. In particular, a prob- 
lem which waits for a thorough illumination is the na- 
ture of multiparticle entanglement . There are several 
aspects of quantum multiparticle correlations that have 
been investigated recently. One example is the investi- 
gation of intrinsic A'-party entanglement (i.e, generaliza- 
tions of the GHZ state jlj]). Another is the realization 
that in contrast to classical correlations, entanglement 
cannot freely be shared among many objects. 

Coffman et al. have recently studied a set of three 
qubits, and have proved that the sum of the entanglement 
(measured in terms of the tangle) between the qubits 1 
and 2 and the qubits 1 and 3 is less than or equal to the 
entanglement between qubit 1 and the rest of the sys- 
tem, i.e. the subsystem 23. Specifically, let us define the 
bi-partite concurrence p6[ of a two-qubit system in the 
state gjk to be: 



5 



c. 



jk 



G{Qjk) ■= max{0, Ai - A2 - A3 - A4} (6.1) 



where the AiS are the square roots of the eigenvalues of 
the matrix R = Qjkic^y ® ^v)i0jk)* (cry (g) ay) Hsted inde- 
creasing order. The tangle is equal to the square of the 
concurrence, i.e. Tjk — [Cjk]'^. Using this definition we 
can express the Coffman-Kundu-Wootters (CKW) 
inequality as 



C12 



Ci3 < 



(23) 



(6.2) 



In the same paper they conjectured that a similar in- 
equality might hold for an arbitrary number, N , of qubits 
prepared in a pure state. That is, one has 



N 

E 

k=Uk^j 



^j.k 



< c 



3,J 



(6.3) 



where the sum in the left-hand-side is taken over all 
qubits except the qubit j, while C^- debotes the con- 
currence between the qubit j and the rest of the system 
(denoted as j). 

Several interesting results in the investigation of var- 
ious bounds on entanglement in multi-partite systems 
have been reported recently. In particular, Wootters jT^ ] 
has considered an infinite collection of qubits arranged 
in an open line, such that every pair of nearest neighbors 
is entangled. In this translationally invariant entangled 
chain the maximum closest-neighbor (bi-partite) entan- 
glement (measured in the concurrence) is bounded by the 
value 1 /\/2 (it is not known whether this bound is achiev- 
able) ||l^. Later Koashi et al. ||l^ considered a finite 
system of N qubits in which each pair out of N{N — l)/2 
possible pairs is entangled (the so-called web of entangle- 
ment). It has been proved that the maximum possible 
bi-partite concurrence in this case is equal to 2/N . Diir 
^ 1 considered other possible inequalities associated with 
variously entangled qubits in multi-partite systems. 

Within the context of our investigation it is very natu- 
ral to ask, what is the nature of the entanglement created 
during the process of homogenization. In this section 
we will address several questions related to this issue. 
Firstly, we will study the bi-partite entanglement be- 
tween the system qubit and the reservoir qubits, and then 
we will analyze entanglement between reservoir qubits 
which is induced by the interaction with the system 
qubit. We will show that CKW bounds are saturated, 
that is the -I- 1 qubit state created by a sequence of 
partial-swap operations in th e ho mogenization process 
satisfies the inequality in Eq. ( x3 ) as an equality. 



A. Bipartite concurrence 



in) 

Let us consider the concurrence Cjj, between the j— th 
and fc— th qubits (irrespective whether these are reservoir 
or system qubits) after the n— th interaction, assuming 



that initially the system was in the state g and the reser- 
voir qubits were in the state ^. Without loss of generality 
we shall always assume that j < k. The value j — de- 
notes the system qubit and j — 1,2, ... ,N denote the 
qubits of the reservoir. 

The reduced density operator g^^^} describing the two 
qubits under consideration is given by the expression 



g%^^Trjj:{Ur....U,[g®e'']) , 



(6.4) 



with Ui[a] — PiaP^ , where Pi is the partial swap opera- 
tion acting between the syst em qubit and the l-ih qubit 
of the reservoir [see Eq.( |2.3| )]. The line over the indices 
in the trace formula denote the partial trace over all sub- 
systems except the ones with the line over them. 

Using the definition (6.1) of the concurrence, it is triv- 
ial to see that C, 



(n) 
jk 



for i, k > n, that is, the qubits 
which have not interacted yet are not entangled. On the 
other hand, a general expression for the concurrence is 
difficult to derive, so we concentrate our attention on a 
special case, when the reservoir is initially in a pure state 
1^) while the system qubit is in an arbitrary state g. 

Following the homogenization scenario the system 
qubit after (fc — 1) interactions is in the state glj' ^\ 
that can be expressed in the basis {|^), l^"*")} as 



(fc-i) 
£•0 = flfc- 



,ie)(ei + (i-afe-i)ie^)(e' 



(6.5) 



After we apply the fc-th partial swap operation between 
the system and the fc-th reservoir qubit we find the bi- 
partite density operator in the matrix form (in the given 
basis) 

(flfc-i cbk-i isbk-i ^ 

6fc_iC (l-afc_i)c2 isc{l-ak-i) 
-isbl_^ -iscil-ak-i) s2(l-afc_i) 
00/ 

(6.6) 

where c — cos rj, s — sin rj. 

The matrix, R, constructed from g'^f} has only one 
nonzero eigenvalue 4c^s^(l — afe_i)^. This implies for 
the concurrence 



d^l^ = 2cs(l - afc„i). 



(6.7) 



From Eq. ( p.6[ ) we find the recurrence formula for the pa- 
rameters ak 

Gk^ak-ic' + s^ = l~c^''{l~ao), (6.8) 

from which we obtain 

C^^^ =2csc2('=-i)(l-ao) (6.9) 

where qq :— {^\g\£,), and Cq^'' is the concurrence mea- 
suring the entanglement between the system qubit and 
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k—th reservoir qubit just after their join interaction (i.e., 
it is supposed that the system qubit has interacted all 
together just k times). We can conclude that the sys- 
tem qubit is entangled with fc-th reservoir qubit. On the 
other hand we can ask whether this entanglement will 
persist after the system interacts later with other reser- 
voir qubits. In order to make the discussion simpler we 
will study a particular case when initially the system is 
in the state 1 1) while the reservoir qubits are in the state 
|0). Nevertheless, prior to this task we study another as- 
pect of multi-partite entanglement within the context of 
homogenization. Specifically, we will study how a given 
qubit is entangled with the rest of the system. 



B. One qubit vs rest of the system 

In the case of pure niulti-qubit states one can define a 
measure of the entanglement between a single qubit and 
the rest of the system with the help of the determi- 
nant of the density operator of the specific qubit under 
consideration. In particular, let us begin the homogeniza- 
tion process with the system and the reservoir qubits ini- 
tially in pure states. After n partial swaps the jth qubit 



is in the state g 



(n) 



Tr-{L(„..M,U)\0^^^^). The de- 



gree of entanglement between the jth qubit and the rest 
of the system is given by the expression 



C. The case |i/')o = |1) and = |0) 

In order to have a deeper insight into the problem of 
entanglement induced by the homogenization process, let 
us consider a specific initial state of the system and the 



reservoir: 



|0) and 10, 



11). In this case we find 



for the tangle between the system and the rest of the 
reservoir qubits after the n-th interaction the expression 



J") 



4deteJJ = 4c^"(l 



(6.13) 



since t^^"^ = (1 - c2")|0)(0| + c2"|l)(l| [cf. Eq. {^]. 
It is clear from this expression that as n — > oo the de- 
gree of entanglement between the system and the reser- 
voir is monotonically decreasing. On the other hand, the 
state of j-th qubit after the interaction with the system 



qubit is 



c2|0)(0| = (l 



-S2c2(j-l))|0)(0|- 



s2c20-i)|l)(l| [cf. Eq. (U)] from which it follows that 



^0") ^4^2^20-1) ^2^20-1)) 



(6.14) 



In other words, after its interaction with the system qubit 
the j-th qubit is constantly entangled with the rest of the 
system. These simple examples illustrate the more gen- 
eral conclusions presented in the previous paragraph. 
Let us turn our attention to the bipartite concurrences 
With the given initial conditions, we easily find the 
state vector describing the whole system after n interac- 
tions: 



M) - 



C, 



:= AAetg) 



(«) 



where is the tangle, which is equal to the square of 
the corresponding concurrence. 

Obviously, for the jth qubit of the reservoir, the tan- 
gle is zero until it interacts with the system qubit. After 
the interaction its value remains constant, irrespective 
of the further evolution of the system qubit during the 
homogenization process. This means that 



.(") 





4det^^ 



if 
if 



n<j<N 
j <n<N 



(6.11) 



In order to justify the last equation we note that all 
measures of entanglement remain unchanged under local 
unitary transformations, and that all transformations Uk 
(except the jth one) are local with respect to the parti- 
tion j © j (where j denotes all qubits except jth) . 

The tangle between the system qubit and the reservoir 
is given by the expression 



(6.10) |*> = c"|l)o ® |0)«^ + V |1)/ ® |0) 



1=1 



We recall that N is the total number of reservoir qubits, 
and that we have assumed that j < k. The state jO)**^" 
denotes all qubits except the qubit I in the state |0). Trac- 
ing over the appropriate subsystems we find the density 
matrices for j < n < k, 



For k < n we find 



Jj) 



(6.16) 



(«) 



and 



„2n 



^2^2(fc-l) 



|00)(00| 



+s2c2('=-i) |01) (01| + s^c^O-i) 1 10) (10| 
+s2c^+'=-2(c + is)fe-J|01)(10| 

+ s2cJ+'=-2(c-is)fc-J|10)(01|, 



(6.17) 



Jn) _ 



in) 



(6.12) 



from which it follows that the shared entanglement be- 
tween the system qubit and the whole reservoir depends 
on the total number of interactions n, unlike in the case 



( |6.11 ) of the reservoir qubits. 



(n) 



J2 s'c2('-i)|00)(00| 
+c2"|10)(10| + s2c2('=-i)|01)(01| 



+isc' 



n+fe-1 



(c + zs)"-^|01)(10| 
(c-2s)"-'^|10)(01| 



(6.18) 
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which determines the values of the concurrences. The In the case n < j all the cj^'' vanish. That means 
corresponding eigenvalues of the matrices R^^\ con- 
structed from the density matrices gpl' (in the case 
n > k), are 



Sj{n) = = rj"\ If n > j then 



eig( 



(i?J^) = {4,sV(-'"+'=-2), 0,0,0} 



eig(i?(",)) = {4.V("+^-i), 0,0,0}. 



(6.19) 



The squ are roots of these eigenvalues are the A^s in 
Eq. (|6.l[). For the concurrences we find 



C 



(") 

jk 



for n < A: < iV 

2s2^j+k-2 fQj. < „ < TV 



in) 



Ok 





2sc" 



for 
for 



n<k<N 
k <n< N 



(6.20) 



(6.21) 



We see that the concurrence between any two qubits of 
the reservoir is zero until both of them have interacted 
with the system qubit. Then the concurrence rises dur- 
ing the interaction to a new value and remains constant 
in the subsequent evolution. On the other hand the value 
of the entanglement between the system qubit and any 
qubit from the reservoir becomes nonzero after their joint 
interaction, but then it tends back to zero. 

This means that the system qubit acts as a "media- 
tor" of entanglement between the reservoir qubits which 
have never interacted directly. It is obvious that later the 
two reservoir qubits interact with the system qubit, the 
smaller the degree of their mutual entanglement. Nev- 
ertheless, this value is constant and does not depend on 
the subsequent evolution of the system qubit (i.e., it does 
not depend on the number of interactions n). 

Once we have derived expressions for the bi part ite con- 
currences, we can verify the CKW inequality ( |6.3[ ), which 
in our notation takes the form 



N 
k=l 



c. 



jk 



< 



c 



in) 



(6.22) 



Firstly, let us consider the entangleme nt of the sys- 
tem qubit with the reservoir. Using Eq. ( |6.21 ) we can 
explicitly evaluat e the expression for the left-hand-side 
of the ineq uality ( 6.22| ) and we can compare it with the 
expression ( 3.13| ) representing the right-hand-side of this 
inequality. We find that: 



^o(n) = E[QrP = 4c^"(l 



(6.23) 



k=l 



which means that the bound on the bi-partite entangle- 
ment between the system and the reservoir qubits is sat- 
urated and the two sides are equal. 

In fact, this property is also valid for the reservoir 
qubits. So, let us consider a qubit j of the reservoir. 



k=l 



■2k 



\^k=l k=j + l 

= 4s2c2("+J-l) + 4s4c2(j-2)(c2 - s^C^^' - C^^^+l)) 

= 4sV(^-i)(l-sV(j"-i)). (6.24) 
In the calculation we used the equality 

J — 1 n n 

fc=l k=j+l k=0 



Comparing this result with Eq. ( |6.14 ) we obtain again 
the equality in Eq. (3.22) 



,2„20-l)l 



for n< N 
for j < n < N 

(6.26) 



To understand in more detail the meaning of above ex- 
pressions, let us consider the entanglement in the limit 
— > oo of a very large number of qubits in the reser- 
voir. We have to be careful with the definition of this 
limit. Let us first recall the definition of homogenization. 
We want to obtain homogenized qubits in states within 
some (^-neighborhood of the reservoir's state (|0)(0| in 
our case). In Section we showed that if we have a 
large number of qubits N , we can achieve a n ar bitrarily 
good homogenization, since in the bound (5.5) we can 
let i5 ^ 0. In turn, the bound (5.2) means that (S — > is 
obtained for s ^ 0. The behaviour of the expression c^^ 
in this limit is as follows: Since s ^ 0, then c 1, but 
still c < 1, th e refor e lim/v->oo c^ ^ = , too. Now, looking 
at Eqs. (|6!20|) , (|6^) , (|6^) and ( |6^ ), we see that in the 
limit A'' 00 all the concurrencies vanish. Therefore, 
the shared entanglement between any pair of qubits is 



zero in this case, i.e. ImiN^oo C, 



(N) 
jk 



0. Also the en- 



tanglement shared between a given qubit and rest of the 
homogenized system, expressed in terms of the function 
Sk{N) is zero: 



N 



lim Sk{N) = 0; 



0,1, 



.A^ 



(6.27) 



So, that's how we define the limit A'^ —^ oo: first, we as- 
sume that the system qubit has interacted with all the A'^ 
qubits in the reservoir, for N finite; then, we let A'^ go to 
infinity, always assuming that we make the best possible 
homogenization according to the bounds of Section 

As a first result, it is instructive to realize that in the 
limit N ^ oo (when (5^0) the functions So{N) and 
Sj{N) are such that 
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lim 



So(N) 
Sk{N) 



= 1, j = l,...N, 



(6.28) 



which means that the entanglement of the system qubit 
with the reservoir is the same as the entanglement of an 
arbitrary reservoir qubit to the rest of the homogenized 
system. This reflects the fact that not only states of 
individual qubits are the same but also the amount of 
entanglement between each of the qubits and the rest of 
the system are equal (see Fig. ||). 



o- 



o 



o 



o 





u, 



01 



U, 



02 



u, 



03 



FIG. 3. In this figure we schematically describe the pro- 
cess of entanglement between the system qubit and the reser- 
voir qubits via the homogenization. The initial state of the 
whole system is shown in the left part of the figure: We have 
the system qubit denoted by the black circle, while the reser- 
voir qubits are denoted as white circles. After the interaction 
(Uoi) between the system and the first reservoir qubits a cor- 
responding change of states (represented in different degrees 
of the grey color) and establishment of the entanglement (rep- 
resented by the thick black line) is exhibited. After the inter- 
action {U02) with the second reservoir qubit, a three particle 
entangled state is created, with various degrees of bi-partite 
entanglement (represented by black and grey lines, where the 
grey line corresponds to the entanglement between reservoir 
qubits which have not interacted directly). In the right sec- 
tion of the figure we see the situation after the interaction 
(f-^03) of the system qubit with the third reservoir qubit. All 
qubits are now entangled, black lines describe the entangle- 
ment between the system and the reservoir qubits which is 
established due to the direct interaction, while grey lines cor- 
respond to the entanglement between reservoir qubits induced 
by the interaction with the system qubit. 

In spite of the fact that the pairwise entanglement be- 
tween qubits in the limit N 00 tends to zero, the 
information about the initial state of the system qubit 
is distributed among the homogenized qubits. Thus we 
have infinitely many infinitely small correlations between 
qubits and it seems that the required information is lost. 
However, as N goes to infinity we have infinitely many 
qubits and the information redistributed among them has 
to be vanishingly small. If we sum up all the mutual con- 
currences between all qubits we obtain a finite value 



lim 



JV 

X 

j<k 



1 



N 



N—^QO Z ^ — ^ 



(6.29) 



This supports our argument that the information about 
the initial state of the system is "hiden" in mutual cor- 
relations between qubits of the homogenized system. In 



the concluding section of the paper we will study how 
this information can be recovered. 



To clarify the meaning of Eq. ( 3.2£ ) , we recall the re- 
cent results of Koashi et al. [|l8|. These authors have 
considered a system of N qubits composing the weh of 
entanglement. That is, each of N{N — 1)/2 possible pairs 
of qubits is entangled, while the degree of entanglement 
is equal for all pairs. It has been shown that the maximal 
degree of pairwise entanglement in the web of entangle- 
ment is given by C*!^' — 2/N, that is the maximum 
tangle is r = i/N^. Given that there are N{N - l)/2 
possible pairs we find that the total value of the pairwise 
tangle is 



lim 

N^oo 



iV(iV- 1) 4 
2 iV2 



(6.30) 



which is the same value as found in the homogenized 
system under consideration. 



VII. CONCLUSIONS AND DISCUSSION: 
APPLICATIONS OF HOMOGENIZATION 

In this paper we have shown that one can choose a uni- 
tary transformation that exchanges information between 
a system qubit and a qubit from a reservoir, which, when 
applied sequentially to the system and each qubit in the 
reservoir, will generate an evolution that has the resevoir 
state as a fixed point. In fact, the state of the system 
qubit and those of the reservoir qubits become the same. 
Moreover, this unitary transformation, which we call the 
partial swap operation, is the only transformation, which 
is independent of the initial states of the system and the 
reservoir qubits, that will accomplish this. 

This result is interesting per se since it allows us to un- 
derstand in greater detail the dynamics of open systems 
ig. It is also a nontrivial fact that the partial swap op- 
eration applied to the system qubit and a set of reservoir 
qubits allows us to realize an arbitrary contractive map 
of the system qubit . 

On the other hand the results presented in the pa- 
per can be used in the context of quantum information 
processing. Specifically, quantum homogenization can be 
utilized for quantum cloning and in a protocol realizing 
a quantum safe with a classical combination. 



A. Quantum cloning 

It is well known that unknown quantum states cannot 
be copied perfectly. Specifically, Wootters and Zurek |2^] 
have presented a very simple proof that a perfect cloning 
transformation for unknown quantum states is impossi- 
ble. The ideal quantum cloning scenario would look as 
follows: The quantum doner is initially prepared in a 
state \S) which does not depend on the unknown state 
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of the input qubit which is going to be cloned. In ad- 
dition, a qubit onto which the information is going to be 
copied is available. This particle is prepared in a known 
state denoted as |0). The perfect copying transformation 
C should be of the form 



\m)\s)^\m)\s') 



(7.1) 



From the linearity of quantum mechanics it follows that 
the cloning transformation (7.1) is not possible. 

Even though ideal cloning, i.e., the transformation 



(7.1), is prohibited by the laws of quantum mechanics 
for an arbitrary (unknown) state \ip), it has been shown 
that it still possible to design quantum doners which op- 
erate reasonably well pl| . These quantum doners have 
been specified by the following conditions: 

(i) The state of the original system and its quantum 
copy at the output of the quantum doner, described by 
density operators p^""*-* and pj""*'' , respectively, are iden- 
tical, i.e., 



^(out) 

Pi 



^(out) 
P2 



(7.2) 



(ii) If no a priori information about the m-state of the 
original system is available, then it is reasonable to re- 
quire that all pure states should be copied equally well. 
One way to implement this assumption is to design a 
quantum copier such that the distances between den- 
sity operators of each system at the output p^""*-* (where 

j = 1, 2) and the ideal density operator pj*'^^ which de- 
scribes the in-state of the original mode are input state 
independent. 

(iii) Finally, it is also required that the copies are as 
close as possible to the ideal output state, which is, of 
course, just the input state. This means that the quan- 
tum copying transformation has to minimize the distance 
between the output state p^°"*'' of the copied qubit and 

the ideal state p'j''^\ 

It has been shown by various authors that quantum 
doners satisfying the above conditions do exist |^ , p2| . 
Recently, experimental realizations of these quantum ma- 
chines have been reported as well [ p3|j2^ ] 

However, this is not the only approach to quantum 
cloning; one can formulate the problem from a slightly 
different perspective using the ideas of quantum homog- 
enization. First, one can lift the condition (7.2) that the 
qubits at the output are in the same state, that is, it 
can be assumed that the qubits at the output are in the 
states which are similar, but not identical. The second 
condition which might be lifted is that the "blank" qubit 
is initially in the known state |0). We can instead assume 
that both the input state of the original and that of the 
"blank" are unknown. If this point of view is adopted 
then the quant um homog eniz ation as characterized by 
the conditions (1.3) and ( |l.4| ) can be successfully used 
for approximate cloning. Specifically, in this scenario the 
reservoir qubits play the role of originals, that is, it is the 



state ^ we want to copy, while the system, which is sup- 
posed to be homogenized, (this system is initially in an 

unknown state g'g '' ) plays the role of the "blank" system 
onto which the information is going to be copied. From 
the description of quantum homogenization we see that 
quantum cloning in this context is a process in which we 
start with N reservoir qubits, all in the same state ^, 
and we end up with -I- 1 qubits in states which are very 
close (how close depends on the value of N) to the state 
^, so we have performed a version of — > A'^ + 1 cloning 
on the reservoir state. 



B. Quantum safe with a classical combination 

After the system qubit is homogenizaed it is in the 
same state as the reservoir qubits, so we can ask: What 
happened to the information encoded in the initial state 
of the system qubit? Is it irreversibly lost? Certainly not, 
because we considered only unitary transformations, and 
that means that the information encoded in the initial 
state of the system qubit is not lost but it is transferred 
into quantum correlations between all of the qubits. The 
parameters characterizing the state of the system are 
transformed into parameters determining the entangle- 
ment shared among the system and reservoir qubits. One 
question is whether the initial state of the system qubit 
can be recovered. 

The process of homogenization is described by a se- 
quence of unitary operations. Consequently, it can be 
reversed: That is the homogenized system can be "un- 
wound" and the original state of the system g^"' and the 
reservoir ^ can be recovered. Perfect unwinding can be 
performed only when the A^-f 1 qubits of the output state 
interact, via the inverse of the original partial-swap op- 
eration, in the "correct" order. The system particle must 
be identified from among the A^ -f 1 output qubits, and it 
and the reservoir particles must interact in the reverse of 
the order in which they originally interacted. Therefore, 
in order to unwind the homogenized system, the classical 
information about the ordering of the particles is vital. 
Obviously, if there are at the output A^-|-l particles, then 
there exists (A^ +1)! permutations of possible orderings, 
only one of which will reverse the original process. The 
probability to choose the system particle, which is in the 

state g^^\ correctly is 1/(A^ +1)- Even when this parti- 
cle is chosen succesfully, then there are still A^! different 
possibilities in choosing the sequence of interactions with 
the reservoir qubits. If one has no knowledge about the 
output particles, the probability of successfully unwind- 
ing the homogenization transformation is 1/(A^+ 1)!. As 
we shall see, if at the beginning of the unwinding process 
the reservoir particle is chosen incorrectly then the whole 
process leads to a completely wrong result. 

Therefore we can consider the quantum homogeniza- 
tion as a process that generates a combination to a quan- 
tum safe. The combination is the sequence in which the 
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reservoir particles interact with the system particles, and 
the object in the safe is the initial state of the system par- 
ticle. The combination consists of classical information, 
and the object in the safe consists of quantum informa- 
tion. The security here is given by the homogeneity of 
the final ensemble; it is difficult to distinguish among 
the particles by measuring them. The unwinding process 
can be performed reliably only when the combination is 
available. An important aspect of this scheme is that if 
one has tried one possible unwinding of the state, and 
measured the result to gain some information about it, 
it is not possible at that point to try to unwind it in a 
different way. That is, the nature of quantum mechanical 
measurement prevents repeated unwinding procedures on 
the same homogenized set of particles. 



60000 
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FIG. 4. The result of the unwinding process with a 
trial-and-error strategy, when the system qubit is correctly 
chosen from a set of 10 homogenized particles. We plot a 
histogram representing number of reconstructed states of the 
system qubit with s falling into the bin with « = z„ ± 0.05. 
There are altogether 9! sequences we have checlted and just 
one results in a correct reversal of the homogenization process. 

To illustrate the above protocol let us assume that we 
begin with the system qubit in the state |1) and nine 
reservoir qubits in the |0). After quantum homogeniza- 
tion we try, randomly, to unwind the process. Let us 
assume that we are lucky and we have chosen the first 
qubit in the unwinding process correctly, that is, it is 
the original system qubit. Even with this good start, 
we have to find the rest of the combination, the proper 
sequence of the reservoir qubits, in order to completely 
"open" the quantum safe. Here we adopt a trial-and- 
error strategy, and we test all possible permutations of 
the reservoir qubits. Obviously, in this case just one se- 
quence is correct, i.e. will result in opening the quantum 
safe and recovering the system state. All 9 ! = 362880 
possible permutations of the reservoir qubits which were 
tested. Since we have chosen the states of the sys- 
tem and the reservoir qubits to be two orthogonal basis 
states of a single qubit, we can parameterize the recon- 
structed system state with just a single parameter s, i.e. 

Punwound = ^(^+^''^^1= ^ |0) (0| + ^ 1 1) (1 1 , such that 
— 1 < z < 1. In Fig. Hwe plot the histogram represent- 



ing the number Ng of reconstructed states of the system 
qubit with s falling into the bin with 2: = z„ ± 0.05. We 
see that a randomly chosen combination will not open the 
quantum safe. In fact, most of the reconstructed states 
are within the interval —1 < z < 0, i.e. between the 
reservoir state and the completely random state. 
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FIG. 5. The result of the unwinding process with a 
trial-and-error strategy, when the system qubit is chosen in- 
correctly from a set of 10 homogenized particles. In the fig- 
ure we represent results of 9 ■ 9! random unwindings. None 
of these sequences lead to the correct reversal of the original 
homogenization. We plot the histogram representing number 
of reconstructed states of the system qubit with s falling into 
the bin with z — Zn i: 0.05. 

Let us now consider what happens when we choose 
the wrong qubit as the system qubit, i.e. what we have 
chosed as the system qubit was, in fact, one of the orig- 
inal reservoir qubits. As can be checked explicitly, in 
this case there is no way to correctly unwind the homog- 
enization process. Obviously, with no prior knowledge, 
the probability to chose an incorrect system qubit from 
a set of -I- 1 homogenized qubits is N time larger than 
the probability to chose the system qubit correctly. In 
addition, there are N ■ N\ different sequences for the un- 
winding procedure in this case and none of them results 
in the initial state. In Fig. ^ we plot the results of these 
unwinding procedures for the same choice of the initial 
states as in the previous case. 

We can conclude that the process of quantum homog- 
enization can be unwound (i.e. reversed) if and only if 
classical information about the order of reservoir qubits 
is available. If this information is discarded, the process 
becomes irreversible even though the overall dynamics 
is unitary. This irreversibility can be used to protect 
quantum information. A detailed analysis of the security 
of the protocol we have proposed to do this remains to 
be done, but the example we have treated numerically 
strongly suggests that quantum information protected in 
this way is very secure. 
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